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1. Introduction 

Lattice gauge theories^'^ have allowed a deep understanding of the contin- 
uum behavior of non-abelian gauge theories. It was not thinkable, before the 
introduction of such a powerful non-perturbative tool, to have such a good 
control of properties like quark confinement and of the gluonic and fermionic 
mass spectrum. 

The string tension, as defined in the quenched theory, which we will mainly 
discuss in the following, is maybe the main non-perturbative quantity. We 
use it to learn that the theory is confined, and to gather non-perturbative 
information about the behavior of quark pairs. 

Monte Carlo simulations have been described in a very complete manner 
at various stages of their developments, both in reprint collections"^, in books 
dedicated to the subject and in reviews®'^'^°. The Lattice conferences 
proceedings provide an invaluable source of information. We will not quote 
them one by one, but address the reader to the yearly reviews on the string 
tension and, more in general, to the detailed contributions contained in the 
proceedings. 

Also thanks to the existence of such reviews here we will not have to 
condense large series of numbers, but we will just try to reconstruct a general 
frame, and to stress a few main ideas, and accomplished and yet to come 
developments. We will try here to give the basics that a beginner would 
need to approach the field, to stress a few ideas that look important to us, 
and to indicate the direction of the most promising developments. Nothing 
more. 

In section (2) we will introduce Wilson loops, and their use to measure the 
string tension. We will discuss there the first string tension measurements. 
In section (3) we introduce Polyakov loops. In section (4) we discuss the 
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continuum limit, and how one gets sure that she is measuring a quantity 
that does not depend anymore on the specific (lattice) discretization. In 
section (5) we discuss the measurements of the quark anti-quark potential. 
Since this is a crucial physical point we will go here in some detail. 

Maybe the main point of the review is in section (6), where we discuss 
a variety of methods to measure the string tension. This has been one of 
the crucial developments which have allowed to use the computers to extract 
non trivial physical predictions. In section (7) we discuss about universality, 
and, again, about which is the best to the continuum limit. In the section 

(8) we discuss, in some sense, about finite size effects. It is interesting that 
the picture of a confinement given by a thin fiux confined to a tube can lead 
to quantitative predictions of the size of the non-dominant terms. In section 

(9) we draw our conclusions. 

2. Wilson loops 

In this section we will discuss the use of Wilson loops to measure the 
string tension in non-abelian gauge theories. This is the simplest framework 
in which non-perturbative quantities can be computed (at least numerically) 
with good precision. Wilson ideas^ find here a practical implementation, and 
become operative tools. In the following we will discuss mainly about the 
quenched theory (where fermions do not appear in the functional measure). 
Here fermions are external sources, which we use as a probe to test the in- 
trinsic confining properties of the gauge field configurations which contribute 
to the path integral. 

One of the main goals of a non-perturbative formulation of Gauge The- 
ories is trying to understand the phenomenon of quark confinement. In a 
non-abelian gauge theory color sources are confined in color-singlets, which 
cannot be separated. At a qualitative level such confinement is understood 
by assuming that the color electric fiuxes emanating from the quarks are 
squeezed into a string-like configuration. Such a gauge field configuration 
has a constant energy density per unit length, i.e. E = ar. This energy den- 
sity (T is called the string tension, and it is the fundamental physical quantity 
for a pure gauge field configuration (in the full theory the string can break 
by creating a quark-anti-quark pair, and we have to use another criterion 
to check confinement). Verifying such a scenario has been one of the major 
tasks of numerical simulations of lattice gauge theories. 

The fundamental quantity one defines when formulating Lattice Gauge 
Theories is the Wilson loop (W^), which will also be one of the basis of our 
string tension measurements. Let us consider a rectangular closed path 7 in 
(4c/ in our case) Euclidean space-time, extending for a length x along one of 
the 3 spatial directions and for a length t along the time axis. For the SU (N) 
gauge theory on defines by (W^) the expectation value of Tr where 
is the transport operator along the path 7 (in the fundamental representation 
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of the gauge group). On a 4 dimensional hypercubical lattice is simply 
given by the product Ui^i^Ui^i^....Uij^i^ of the (oriented) dynamical variables 
associated with the links joining the consecutive lattice sites one finds along 
the paths. 

In the limit of large separation (W^) describes the change in the vacuum to 
vacuum transition amplitude induced by the presence of an external current. 
The current corresponds to the creation at some initial time of a pair of static 
sources at separation x, the propagation of the sources for a time interval 
t and the final annihilation of the sources. Wilson loops are a probe for 
quarks. From their large separation behavior one learns about the behavior 
of hypothetical (infinitely heavy) quark when put in an typical equilibrium 
gauge field configuration. As we will see a confining scenario implies an area 
law decay of the Wilson loops expectation values. 

If we expand over a complete set of eigenstates of the Hamiltonian (in 
presence of an external source) (W^) can be expressed as: 

{W,) =Y.Cn e-^"* , (2.1) 

n 

where the C„ are the matrix elements for the creation of the sources from the 
vacuum and the exponential factors are real because of the Wick rotation to 
imaginary time (see for example ^ and ^^). 

For large t the state with the lowest energy eigenvalues Eq dominates the 
statistical sum. From the large t and x decay of (W^) one can infer the 
behavior of the energy density* 

Typically one finds that at large x and t 

{W^) ^ ^-^^t-m(^+t)+... ^ ^2.2) 

with an area and a perimeter term. The perimeter contribution e"'"^^"'"*) is 
suggested from the form of the self-energy term in the perturbative compu- 
tation of the heavy quark potential (see for example ^). Now one measures 
numerically the expectation value of W^, by the standard sampling tech- 
niques. If asymptotically one finds that the Wilson loop decays according to 
an area law, i.e. 

{W^) ^ e"""* , (2.3) 

*It should be always clear that the definition of the string tension obtained from the 
asymptotic behavior of the quark-quark potential is in principle different from the one 
obtained by using directly the energy of the electric fiux tube, and that the two have to be 
equal only asymptotically. For example the size of the string fiuctuations can be different 
in the two cases. See for example ref. ^'^ for an interesting discussion of the subject, and 
references therein. 
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than the theory is confined. The effective potential is asymptotically linearly 
increasing, and one estimates the value of a by such a decay rate. On the 
contrary a theory with free particles gives raise to a perimeter decay law. 
That is for example what happens for electrons in the U{1) gauge theory (in 
the continuum limit). 

Let our rectangular loop 

{W,) = W{I,J) (2.4) 

(on a hypercubical lattice in four Euclidean space-time dimensions) have a 
I X J size [x = la, t = J a, where a is the value of the lattice spacing). As 
x — > oo W{I , J) can be expressed in terms of the potential V{x) ^ ax = ala 
as 

W{I,J)^e-''^'^-' . (2.5) 

A lattice calculation of W{I , J) will not produce a directly but a dimension- 
less (since we always measure pure numbers) function A", depending on the 
lattice coupling constant g through the lattice spacing a, related to a by 

K{g) = (To^ . (2.6) 

In the following we will also use the variable fi = 2N j g"^ for the SU (N) gauge 
groups. Section (4) is devoted to study this dependence. 
To measure a Creutz introduced the ratio^"* 

W{I,J)W{I-1,J-1) 

w{i,j -i)w{i ■ ^ ■ ^ 

This ratio has the advantage of making more clear the area dependence of the 
Wilson loops. If W behaves as we have supposed in (2.2), then the Creutz 
ratio has the value 

exp —aa^ . (2-8) 

Obviously for small loops short distance corrections will spoil the behavior 
(2.8). One will have to measure loops that are large enough (in some sense we 
will discuss better in the following) and check a posteriori, self-consistently, 
the presence of a confining behavior and of an area law decay. Technically 
one takes the logarithm of the ratio's (2.7), and defines 

X[l,J)- ^''w{I,J-l)W{I-l,J)- 

When measuring these quantities difficulties arise from the finite extent of 
the lattices. Moreover when W{I, J) becomes too small the error due to sta- 
tistical fiuctuations makes the measurements meaningless: on a finite lattice. 
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for finite / and J, one always estimate a value of x larger then the true value 
of cra^. To extract the string tension one has to use the envelope of the curves 
J)- Let us stress that in the x ratios we have canceled only the dominant 

perimeter behavior of (2.2). Contributions like for example a term e~ (ex- 
pected from the fluctuations of the surface of the path^) do survive, and only 
cancel in the asymptotic limit. We will give a more accurate analysis of the 
quark-antiquark potential behavior in section (5). For a pedagogical review, 
containing explicit examples about the quark potential and the asymptotic 
freedom in continuum quantum cromodynamics, see^. 

First numerical results for the SU (2) and SU (3) gauge theory were pre- 
sented in 14, 15, 16, 17,18, i9_ Computations for different lattice sizes have been 
performed and the dependence of the results over the Fuclidean time exten- 
sion of the lattice has been analyzed. 

From such first string tension measurements it became clear the crucial 
role of the control of finite volume and finite Fuclidean time effects, of the 
contamination from excited states and the distortion of the spectrum due 
to the finiteness of the lattice. In particular we will see that the difficult 
problem of extracting the rate of an exponential behavior out of numbers 
of order one has been overcome not only with more powerful computers, 
but (and maybe this has been the crucial step) also by developing effective 
numerical algorithms. 

3. Polyakov loops 

Polyakov lines, which wrap the periodic lattice in one direction, have been 
first used in the context of Lattice Gauge Theories, in ref. ^° and to mon- 
itor a finite temperature phase transition. Polyakov lines turn out to be, as 
we will discuss in the following, very useful tools to get reliable measurements 
of the string tension. 

After the first string tension measurements, which we have discussed in 
the previous section, larger scale measurements of the Wilson loops, at lower 
values of the inverse coupling constant /3, suggested that the asymptotic 
string tension could assume a value substantially smaller than previously 
determined. 

It was remarked in^'' that a potentially effective way of measuring the 
string tension was based on studying the large separation behavior of corre- 
lation functions of two Polyakov loops (see fig. I). 

Denoting by P[x^y^z) the trace of the product of the link gauge fields 
in the t direction, closed thanks to the periodic boundary conditions (the 
so called Polyakov loops) ^ let us consider in a Ad lattice, with Nf sites in 
the t direction, the correlation functions of two Polyakov loops wrapping 
the lattice in the t direction, separated in the z direction by z sites. In the 
scaling limit, for z large enough and — > oo, one recovers the quark- anti- 
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Figure 1: (a) Generic correlation of Wilson loops of length R separated by 
the time distance T; (b) Correlation of Polyakov loops encircling the periodic 
time direction. 

quark potential by 

(P(0, 0, 0)P(0, 0, z)), -AT.^oo e-"^' ""^'^ . (3.1) 

If in the large z region V{z) increases linearly with the distance, i.e. V{z) 
Kz^ then 

(P(0, 0, 0)P(0, 0, z)), -AT.^oo e-"^' '''^ . (3.2) 

In other words [Nf K) is the mass of the correlation function of two Polyakov 
loops. 

The coefficient {Kz) which governs the decay of the correlation function 
can be seen as the energy Eo[z) of a string of length z which propagates for 
an (Euclidean) time Nf. As for a glueball state (see later), the mass at rest 
of the system can be extracted by summing over spatial planes to project out 
the zero momentum component of the state. Since the straight Polyakov loop 
is already symmetric in the z direction, only an average over the remaining 
two directions is needed. The quantities considered when measuring the 
SU (3) string tension^^'^"^'^'*'^^'^^'^'''^® are correlations of zero momentum sums 
of Polyakov loops: 

C{z) = (^[P(x, y, 0)] E[P(u;, V, z)]), ^iv,^oo,.^oo e-"^' '''^ • (3.3) 

For a lattice of sufficiently small size in the z direction there will be a 
deconfining phase transition (a first order transition for S'[/(3)^^'"^°'"^^'"^^ and 
a second order transition for S'[/(2)^°'^^'"^^'"^"^'"^^). In this gluon plasma phase 
the correlation between two Polyakov loops is given by an inverse power law 
rather than by an exponential because of the zero mass deconfined gluon 
modes. We use fig. 2 to illustrate the situation. 
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Figure 2: The expected form of the energy Eo[z) as function of z: Eo[z) = 
for z < Np where Np is the finite temperature transition in a Np X oo"^ lattice. 
The extracted value of Eo[z)/z = Keff leads to Keff ^ K where K is the 
asymptotic slope of Eo[z) 

The string tension extracted from a finite lattice will be smaller than the 
true string tension given by the asymptotic slope of Eo[z) in z. The presence 
of such finite temperature effect makes clear the fact that Eo[z) cannot be 
linear with z. Considering contributions of different nature is needed for 
achieving a complete and satisfactory description. 

The general requirement needed for obtaining a fair measurement of the 
string tension is to use a lattice large enough so that the contribution from the 
fake quark loops wrapping around the lattice are sufficiently small. This effect 
is connected to the finite temperature phase transition which is triggered 
when such a contribution becomes significant. In fig. 3, from ref.''^, we show 
the minimum lattice function of /3, needed to avoid finite size effects 

for SU{3). 

A complete review of the measurements of SU(2) and SU(3) string tension 
until 1987 can be found in ref.''^; more recent data can be found for SU(3) in 
60,58,69^ for SU(2) in j.ef.i™'^®'63'^9 -^^ ^j-^g references discussed in more 

detail in the following sections. 

An interesting comparison between the string tension measured from Wil- 
son loops and Polyakov loops can be found in references"^'*'"^^. In particular 
ref."^^ shows that using the results of a simplified string model for the confined 
pure gauge field theory it is possible to explain some of the systematic effects 
in the measurement of the lattice QCD string tension. In fact, the finite 
lattice used in computer simulations causes the string tension to be over- 
estimated when using Wilson loops and to be underestimated when using 
Polyakov loops. 
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Figure 3: Minimum lattice size estimated from the critical coupling of the 
finite temperature phase transition as a function of (3 from ref/^. 

Let us finally remark that Polyakov loops method can be generalized to 
the case of multi-quark systems to extract the static energy of the system 
containing quarks and antiquarks"^^'^^'"^''. 

4. Continuum limit, scaling regime and /3-function 

We want now to discuss how from a lattice theory one can recover a contin- 
uum theory. We are interested in continuum QCD, the non-abelian theory 
which we expect to describe strong interactions, and the lattice theory is for 
us a very useful tool, since it provides a gauge invariant non-perturbative 
renormalization scheme for the continuum gauge theory. The lattice theory 
is gauge invariant for all values of the lattice spacing, this is the big bonus 
we get thanks to the Wilson formulation. The fact that we are dealing with 
an asymptotically free theory (i.e. that quarks behave as free at very short 
distances)"^®'"^^ implies that in the lattice model we need to send (/ — > to 
obtain the continuum limit (of zero lattice spacing, a — > 0). The point that 
explains this fact can be visualized by thinking at a finite p/ij/s«ca/ interaction 
range. Continuum limit means that the lattice grid is irrelevant, i.e. in a 
constant physical distance (given by the interaction range) we have more and 
more grid points. But that happens when the interaction becomes weaker, 
i.e. in our theory in the short distance limit. 

In some sense QCD is a happy and easy case, since the critical exponents 
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can be computed in perturbation theory; around g* = 0, the infrared fixed 
point of the theory, perturbation theory and renormahzation group are used 
to understand the critical behavior. Numerical simulations allow to extract 
the non-perturbative content of the theory. Then an analysis is performed 
on the scaling region of the theory, where the perturbative scaling regime is 
expected to hold. 

To give a hint about the behavior of the mass gap of the theory, we use 
the lattice regularization noting that what we are saying would be valid in all 
schemes (but only the lattice scheme is non-perturbative and gauge-invariant 
and will allow the next, crucial step, numerical simulations). 

One requires that a physical quantity with dimensions mass M be a renor- 
mahzation group invariant, i.e not to depend on the lattice spacing in the 
scaling regime which characterizes the approach to the continuum limit: 

dM , , 

= , (4.1) 

where a is the cutoff. Since M has the dimensions of a mass, in the continuum 
limit, where the correlation length diverges in lattice units and there are no 
scales left (the lattice grid disappears in this limit) we have 

M = a-^ig) , (4.2) 

where (/> is a function of the coupling g. Substituting M in (4.1) we see that 
(f)' behaves as 

^ -llfl . (4.3, 

Here fi{g) stands for the well known Callan-Symanzik fi function'*^ 

m = -I . (4.4) 

For non-abelian quantized gauge theories it is possible to compute in per- 
turbation theory the functional form of (j){g) in eq. 4.3 in the vicinity of the 
fixed point g* = 0. This leads, for a pure SU{N) theory, to the result s^'^Q. 

= (3og^ + (3^g' + .... (4.5) 

Oa 

with 

- y (j^), A - -jij^.r . (4.6) 



Using eq. 4.5 the formal integral for (j){g) following from eq. 4.3 i 



is: 
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m = ifJog')^'o^ exp{--fJog'}il + Oig')) . (4.7) 

Then going to the continuum hmit requires the knowledge of a definite rela- 
tionship between g and the lattice spacing a which can be expressed in the 
general form as: 

a = j(j){g) , (4.8) 

where </>((/), as we have seen, is a definite function containing the information 
about the correct scaling behavior of a when g ^ g* , given by eq. 4.7, and 
A is a scale parameter. 

In principle all physical quantities can be expressed in terms of A. One 
goes to the continuum limit by demanding that one (generic) observable of 
dimension — o?, let us say gi, remains strictly constant throughout the process 
of renormalization, i.e.: 

a^a{g) = [q,]'"U.gf'\ (4-9) 

where 4'i{g) measures the observable qi in units of the lattice spacing. The 
scaling properties of the critical point involve that all the function (t)i{g) cor- 
responding to the generic physical quantity g^, of dimension c/i, must behave 
as: 



Ug) ~ c,[ct>{gt , (4.10) 

for g ^ g* , the q being definite constants. Then the continuum values of 
the observables qi are given by: 

q^ = c,A'' , (4.11) 

i.e. all physical observables are expressed in terms of the scale parameter A. 

Then any observable might be used to establish the scale of masses and, 
eliminating A , all other observables are expressed in terms of that one. The 
string tension has traditionally this role. 

A does not have a direct physical significance and in general it will de- 
pend on the scheme of renormalization and on the specific choice of the lattice 
gauge action, within the class of actions leading to the same continuum the- 
ory. A will change if the renormalization scheme is modified. Then in order 
to relate the mainly non-perturbative lattice world to results obtained by 
using perturbative methods (and even, for example, different lattice results 
obtained by using different lattice actions) the lattice scale factor Al must 
be compared with the scale factors defined in the conventional scheme of per- 
turbative renormalization as Ajg5, Amom and so on. Computations of ratio 
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between lattice and perturbative scale parameters can be found in 



42,43,44,45,46 



Let us summarize. Near the critical coupling g* = a physical observable 



F having a dimension of {massY has to behave according to 

1 



F~cA/,Ai = -</)((/) . (4.12) 



a 



where c is a constant which is not computable in perturbation theory and 
(j){g) is given in eq. 4.7 for a pure SU{N) theory. Then the expression for 
Al in SU{3) is: 

It is a very crucial step in a numerical study to examine whether the mea- 
sured physical quantity obeys this scaling law. This is all the point. In 
a numerical simulation we will study the model for different values of the 
coupling constant. After checking that finite size and time effects are under 
control we will have to establish that indeed the measured observable follows, 
as a function of the coupling constant, the scaling behavior. If this is true 
the measured numerical constant c has a physical meaning and represents the 
value of the physical observable we are looking at, in the continuum limit. 

Then the strong coupling results valid for large values of g must be extrap- 
olated to a value of g small enough that scaling toward the continuum limit is 
seen to take place. Also, the extrapolation should provide a self-consistency 
check for the theory, giving indications that no intervening critical points 
make the properties of the system totally different in the strong coupling do- 
main and in scaling domain. Only in this case features such as confinement of 
quarks, which can easily demonstrated in the strong coupling domain, would 
survive the passage to the continuum limit. 

The string tension has played a role of primary importance in the pioneer- 
ing investigations of the scaling of lattice quantities and to verify asymptotic 
freedom. This quantity is regarded as the most fundamental physical quan- 
tity of lattice gauge theory and it has been continuously checked (even as 
a check of the computer program!) by workers in the field when they start 
their simulations. In particular the confining properties of the theory mani- 
fest themselves in non- vanishing string tension a. 

According to the general principles since a has a dimension (mass) we 
expect it will scale as ~ Al^ on the lattice close to the weak coupling limit 
(/ — > 0. The aim of numerical simulation is, therefore, to extract a non 
perturbative quantity: 



c = Al/v^ , 



(4.14) 
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Figure 4: The quantities /) for SU{2) in (a) and for S'[/(3) in (b), as a 
function of l/g^. The envelope of the curves describes the string tension as 
a function of the couphng. 



which connects the scale parameter with the string tension that controls the 
hadron spectrum. The string tension is a concept defined rigorously only in 
the absence of dynamical quarks. In their presence, it becomes energetically 
more favorable for large r to split the string in the middle by creating a 
quark-antiquark pair out of vacuum. 

The first investigations of the string tension were devoted to verify the 
scaling behavior of (t((/) as (/ — > and to use such a scaling behavior to find 
the relation between string tension and lattice scale parameter. 

The studies in the SU (3) Hamiltonian formulation'*''''*® and in SU (2) and 
SU (3) Euclidean formulation'*^'^° are based on the strong coupling expansion 
and relay on some method of extrapolation to probe the behavior of the 
theory as g becomes small. 

At the other hand the first numerical results presented 
based on Monte Carlo simulations ad are independent of any extrapolation 
technique. The evidence for the scaling derived from this data appeared 
particularly impressive (we give in fig. the world famous Creutz figure, so 
suggestive at the time). 

All this first studies are concerned with verifying the scaling behavior of 
(t((/) as (/ — > and with using such a scaling behavior to find the relation 
between the string tension and the lattice scaling parameter. 
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The string tension, defined as the coefficient a of the linearly rising part 
of the potential for large separations of the quark-antiquark pair, in a nu- 
merical calculation can only be extracted from relatively small qq separations 
(in lattice units). Hence if the scale at which this force is seen in nature is 
determined by the size of hadronic matter one must choose the lattice to be 
sufficiently coarse (i.e. the coupling g sufficiently large) for hadron, repre- 
sented by the qq system, to fit in it. But if the lattice is made too coarse one 
cannot expect to extract continuum physics, and the dimensionless string 
cra^ will not exhibit the behavior predicted by the renormalization group. 
Then if one chooses too small a lattice spacing (or coupling constant) then 
one cannot expect to see the asymptotic form of the potential on a lattice 
of small extent; if the lattice spacing is too large one cannot expect to see 
continuum physics. Thus at best one can measure the physical string tension 
in a narrow range of the coupling constant, picking up in this way a narrow 
scaling window. We will see at the end of this review that nowadays the 
situation is far more convenient, and measurements are very reliable. 

We have seen that the /3-function describes the relation between the bare 
coupling constant g and the value of the cut-off, and it has a well defined 
meaning in the vicinity of the fixed point g* = 0. As asymptotic scaling is 
approached the /3-function helps in connecting numerical studies with pertur- 
bation theory. The /3-function describes the way the bare coupling should be 
tuned in order to keep all the physical predictions independent of the cut-off 
in the continuum limit. 

It should be noticed that only the two leading terms of the /3-function 
(3o and (3i are universal. For large cut-off (small bare coupling values) these 
terms dominate and define a universal scaling behavior, the asymptotic scal- 
ing. Outside this region, but still in the continuum limit, the scaling behavior 
is described by the full, and in general, completely unknown /3-function.In 
an SU{N) gauge theory the leading corrections to this behavior are expo- 
nentially small in the inverse bare coupling constant g~^. On the other hand 
the higher order corrections to the universal part of the /3-function are power 
like and are not necessarily small in the region where the cut-off dependent 
corrections are already negligible. There might be sizable contributions to 
the /3-function from non-perturbative phenomena also. Then the non per- 
turbative /3-function tells us how the lattice spacing goes to zero as (/ — > 0. 

On the lattice all dimensionfuU quantities like masses are measured in 
units of the lattice spacing a. In order to take the continuum limit we need 
to know how a scales. One option is to use the 2-loop perturbative result 
provided it is demonstrated that this is valid at values of g where the cal- 
culation (for example of string tension, glueball, etc.) are done; the other 
is to measure the non-perturbative /3-function which quantitative structure 
is necessary to now since the value at which asymptotic scaling sets in is 
not a priori known. One should confirm also that the /3-function approaches 
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its the asymptotic form without undergoing a phase transition, providing 
a continuum hmit with the expected properties of asymptotic freedom and 
confinement. 



5. The q — q potential and improved coupling constants 

When attempting to construct an effective potential to describe the interac- 
tion of quarks one is restricting his focus to heavy fermions. It is only when 
the mass m is large that one can formulate the bound state of qq system as a 
non-relativistic problem, with binding energies computable from a potential 
via the Schrodinger equation. 

General theoretical arguments give the behavior of the spin independent 
potential in the two extreme regions. At large separations (r — > oo), confine- 
ment dominates and the physical picture is the one of a chromo- electric flux 
tube. Here the potential V{r) behaves as a linear function of the distance, 
V{r) — > err. At the other end of the distance scale, i.e. at short distance, for 
small separations of the q — g-pair, the effective potential is expected to ap- 
proach the one gluon exchange result (Coulomb potential). In the continuum 
the perturbative form of the potential (we are dealing with an asymptotically 
free theory) in terms of the short-distance q — g-force A^ parameter^^, known 
by a two-loop computation, is: 

V{r)^-\^, (5.1) 

where the effective coupling Ciqq{r) behaves like for small r. One has 
that 

MO = ^[/3oln(rA,)-' + /3i//3olnln(rA,)-T' , (5-2) 

where the usual coefficients in the perturbative expression for the /3-function 
are given in section (4). The simplest form, in the absence of pair production 
processes, for the full effective potential is 

y(r) = -^^ + ar, (5.3) 
r 

where the two quantities o,qq{r) and a physically represent the scales at which 
an individual term begins to dominate and can be related to the charmonium 
or bottonium spectrum. The study of the potential can be done breaking up 
the separation r between qq into three regions: 1) confining, characterized by 
the linear term at large r; 2) perturbative characterized by a running coupling 
constant; 3) intermediate r values where only phenomenological forms can 
be tried and investigated. 
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The form of the ground state potential, with a Coulomb short-distance 
contribution and a long-range linear piece is well known and helps to account 
for the relative success of potential models in explaining the observed spec- 
trum of heavy-quark mesonic states. These quark-model states do not exploit 
the explicit gluonic degrees of freedom present within the QCD Lagrangian, 
but it is possible to explore this gluonic sector through the lattice gauge the- 
ory approach in all three region of r. In particular the average value W of an 
R X T rectangular Wilson loop is related to the spin independent potential 
V{R) between static color sources on a lattice of spacing a [V{Ra) = V{r)) 
by: 

V{R) = - l\m ^lnW{R,T) . (5.4) 

Thus in principle the knowledge of the potential V{r) allows to determine 
the dimensionless ratio cr/A^ which relates the perturbative scale A to a non- 
perturbative observable such as the string tension a which is, as we have 
seen, the coefficient of the linearly rising piece of the potential. For a proper 
determination of the string tension one should include in the fit the behavior 
of the potential at short distances. If the crossover from the perturbative 
Coulomb potential to the non-perturbative long distance 1/R term is very 
complex a larger than expected intermediate distance Coulomb term could 
be contaminating the values of the string tension. 

Nevertheless the lattice calculations of the potential (where dynamical qu- 
arks are ignored) do not predict a precise functional form: one has to make 
a trial ansatz and use the data from these calculations to fix the unknown 
parameters. Then from a technical point of view the measure of the string 
tension is somewhat problematic, because in practice it is determined by 
extrapolating the heavy quark potential V{r) from a distance r less than 1 
Fermi to large distances. Depending on which analytical form is assumed 
for the extrapolation, one can obtain quite different results. Then the string 
tension is affected by a systematic error which is not easy to control. 

The measurements on the lattice of the string tension and of the potential, 
both related to Wilson loop measures, have parallel histories initiated for the 
potential in53.54,55,56,57_ 

Here we focus our attention on the state of the art measurements of the po- 
tential for SU{3f^'^^'^°'^^ and for SU{2f^'^'^'^^ We will stress the particular 
relevance of the key issue in taking the continuum limit of lattice QCD: the 
question whether the asymptotic scaling behavior holds for physical quanti- 
ties as ^/a/AMs, where^^ ^ 28.81Al for SU{3) and Ams = 19.82Al for 
SU{2). 

A large body of evidence exists that the asymptotic scaling in the the 2 and 
3 color QCD lattice gauge theory is not satisfied at the values of (3 accessible 
at present (up to (3 = 2.85 for SU{2f^ and (3 = 6.8 for SU{3f°. The 
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deviation from the asymptotic scaling behavior is ascribed to finite lattice 
spacing effects. It is possible, however, that at least some part of the deviation 
is due to the presence of perturbative terms of 0{g^) in the ratio \^/^ms 
which can be removed by a redefinition of the coupling constant. Parisi has 
first observed^^ (even before numerical results were obtained) that the bare 
lattice coupling constant = 2N/ (3 receives a large renormalization from 
gluon tadpole contributions, which makes it a poor choice as an expansion 
parameter. Starting from a different approach the authors of ref.^'''^® have 
reached similar conclusions. All together we have a systematic and quite 
crucial improvement of the understanding of the numerical results and expect 
that the scaling behavior of physical quantities will be improved with the use 
of different improved (if so we can say) coupling constant. 

The choice done in ref.^^ is based on a mean-field version of the relation 
between the bare lattice coupling constant g and a renormalized coupling 
constants such as gMs- The expansion parameter gE suggested by Parisi is 
defined introducing in the average plaquette weak coupling expansion: 



an effective coupling in terms of Monte Carlo generated plaquette expectation 
value: 



for which the first-order expansion is supposed to be exact. Further examples 
of different coupling constant schemes can be found in and references 
therein. 

Favourable tests in a number of examples^'''^® show a better agreement 
between perturbation theory and Monte Carlo results which tend to be sys- 
tematically smaller than the perturbative results suggesting that a coupling 
defined from any of the Monte Carlo calculated quantities would yield im- 
proved predictions for the others. 

In this frame the q — q potential is a natural candidate to define improved 
coupling constants. A phenomenological method for estimating the contin- 
uum coupling constant defined by the potential using short distance data has 
been proposed for SU(2) in ref.^"* and extended for SU(3) in ref.^"^'^°. The 
basic idea is to calculate the coupling: 



oo 




(5.5) 



n=l 





(5.6) 



(5.7) 



through numerical simulation and then to convert to the MS renormalization 
scheme warning that the conversion calculated in the one-loop approximation 
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is only applicable when the coupling is sufficiently small, i.e. deep in the 
perturbative short distance regime. 

In ref.^° two approaches are presented for the determination of the QCD 
scale parameter A. First the cutoff parameter Al is calculated from the 
two-loop-expansion equation, following eq. (4.13) 

Then the authors determine the running coupling aqq[R) starting from its 
symmetric discretization in terms of the lattice potential measurements: 



3 „ „ V{R,)-V {R 
~A^'^' R,-R 

whit 



agg{R) = -R,R2^^ , (5.8) 



R = i^l±M (5.9) 



following the ansatz of ref. for the potential: 



V{K) = Vo + (jR-e/R + f/RK (5.10) 

The parameters Vo,(T ,e and / come from a best fit. Then analyzing the 
Q,qq{R) data in terms of the continuum large-momentum expectation for the 
SU (3) running coupling : 

Mi?) = -^[/3oln(rA«)-2 + /3i//3olnln(rA«)-2]-i . (5.11) 

and using the relations A^^ = 30.19A/, and Kj^g = 28. 81 A/, the results are 
converted into the continuum renormalization [MS) scheme. The two meth- 
ods yield consistent values for the Ajy^5-parameter 



Am5 = 0.555i°:°l? , (5.12) 

which is in agreement with the value extrapolated from the scaling of the 
string tension^'''^®. 

This result is substantiated improving on scaling violations by replacing 
the bare coupling with suitable "effective" couplings: the jS^^"^ coupling^^ 
defined by truncating the weak order expansion of the plaquette (measured 

in Monte Carlo simulations) after the first order term, and the scheme 
defined by truncating the same expansion after the second order term. It is 
demonstrated that scaling violations on the string tension can be considerably 
reduced allowing for a safer extrapolation of A/,, as a function of the different 
couplings, to its continuum value; the effective schemes help to decrease 
the uncertainty of this limit and teach that linear extrapolations can be 
misleading: see table 1 from first of references cited in 
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ay/(T 






y^/Af 


5.7 


0.4099 (24) 


124.7 (0.7) 


63.3 (0.4) 


55.7 (0.3) 


5.8 


0.3302 (30) 


112.4 (1.0) 


63.0 (0.6) 


55.6 (0.5) 


5.9 


0.2702 (37) 


102.9 (1.4) 


61.2 (0.8) 


54.3 (0.7) 


6.0 


0.2265 (55) 


96.5 (2.3) 


60.0 (1.5) 


53.4 (1.3) 


6.2 


0.1619 (19) 


86.4 (1.0) 


56.9 (0.7) 


50.8 (0.6) 


6.4 


0.1215 (12) 


81.3 (0.8) 


55.7 (0.5) 


50.0 (0.5) 


6.8 


0.0730 (12) 


76.9 (1.3) 


55.7 (0.9) 


50.4 (0.8) 


oo Lin. 





63.6 (2.4) 


53.1(1.6) 


48.3 (1.4) 


Log. 





54tf5 


ro + 2.6 
'JO .Zj _'j 2 





Table 1: Al as obtained by inserting the bare lattice coupling and A^ ' by 
inserting the [3^^' effective couphngs. Linear extrapolation to a = results 
are displayed in the second to last row and logarithmic extrapolation results 
in the last row 

One observes that the values for A^^ extracted from the running of the 
coupling^'*'^"^'^° tends to increase with the momentum cutoff. The approach 
towards asymptotic limit for different effective coupling schemes are com- 
pared^° in Fig. 5. All that makes clear that lattice perturbation theory works 
better if an improved coupling constant is used. 

In order to overcome the basic difficult encountered in these calculations 
to reach momenta q greater than a few GeV it has been observed™ that the 
running coupling at high energies may in principle be computed once the 
parameters of the theory are fixed at low energy and then is possible to make 
contact with the scaling region where the running coupling is determined by 
the perturbative renormalization group. In this frame it has been proposed 
to use a coupling constant 

«(9) = ^, g=l/i^, (5.13) 

running with the lattice size L. The strength of the gluon interaction at 
high energies can be determined by studying the scaling behavior of such 
an object in small and intermediate volumes. This definition is based on 
a recursive finite-size technique based on renormalization group arguments 
and is composed by two steps. The first is to match lattices with the same 
renormalized coupling but different lattice spacings; the second is to study 
of the evolution of the coupling by changing L at fixed lattice spacing a. 
By combining these steps it is possible to go to the continuum limit and to 
follow the behavior of the renormalized constant g{L) over a large range of 
box sizes L given in unit of the string tension a. 
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Figure 5: Violations of asymptotic scaling for different effective coupling 
schemes. The solid line corresponds to = 5f .9}gA/,. 

Nevertheless the couplings defined through Wilson loop expectation val- 
ues are not a perfect choice since large loops are difficult to compute nu- 
merically (the signal-to-noise ratio is exponentially decreasing with the size 
of the loops) and the perturbation expansion of these quantities to higher 
orders of coupling constant may require an unacceptable amount of work. 
Then g{L) is better defined through the response of the system to a constant 
color-electric background field, or directly by using Polyakov loops''"^''''*. 

This non-perturbative finite-size scaling technique has been used to study 
the evolution of the running coupling both in SU (2)^^''''* and in SU (3)^^ gauge 
theories: one finds again that the perturbative scaling regime and the low- 
energy domain of the theories are smoothly connected with "no complicated" 
transition region. In particular the following relation holds within errors: 

'2aMs{q)\su(2) = ^aMs{q)\su(3) , (5.14) 

where a]^g[q) is determined at momenta q up to lAGeV with an estimated 
precision of a few percent; this relation is consistent with the evolution of the 
coupling up to 3-loop order of perturbation theory and suggests that Naj^g 
is only weakly dependent on the number N of colours. 

Some discrepancies between different numerical measurements of the run- 
ning coupling can be due to non negligible non-perturbative corrections to 
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the matching relations^^ between the different schemes, since the distances 
at which the comparison can be made are rather iarge in physicai units. 
Cutoff effects are hardiy Icept under controi when aqq[R) is determined from 
the heavy quark; potentiai at distances as iow as 2 or 3 iattice distances. A 
summary of various recent analysis for SU (3) contains results all in the 
range: 



to be compared with / Al = 80.0 it 1.4 at /3 = 6.5 extracted by using the 
bare couphng. For SU{2Y^'^^: 



to be compared with / Al = 44.1 it .6 at /3 = 2.85 extracted by using the 
bare coupling. These values are sufficiently far below the values extracted 
from bare coupling to imply that asymptotic scaling to two-loop perturba- 
tion theory is not just around the corner but will only be satisfied with good 
precision at larger /3-values than those currently accessible to lattice simula- 
tion. 

Nevertheless we note that the lattice methods to determine the running 
coupling give an accuracy comparable with that of experimental determi- 
nations of the same quantity for modest energy scales. In particular the 
experimental determination of the running coupling constant of QCD has 
reached a reasonable degree of accuracy''^. Further improvements of lattice 
techniques are of great interest after the discovery that the running coupling 
is well described by the two loop-formula down to a scale of 1 — 2 GeV till to 
be able to predict experimental numbers like a[Mz) or as explained 

in ref.*^^. 
6. Methods 

This section will be devoted to a brief illustration of different methods which 
have favored big progresses towards the study of the continuum limit and of 
large volume lattices in the numerical computation of the string tension. 

Different methods of investigation of lattice gauge theories have been used 
both to measure the glueball spectrum and the string tension: in fact, as we 
have seen in section (3), Nf K is the mass of the correlation function of two 
Polyakov loops and the string tension is related to the ground state of the 
interquark potential. 

We can determine the mass of any state in an euclidean lattice calculation 
by measuring the exponential fall-off of correlation functions of appropriate 
interpolating field operators expressed by the sum over all states that couple 
to them: 




(5.15) 




(5.16) 
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{G{t)G{0))c ^Y.^'n exp [-mj] , (6.1) 

n=0 

that, for t — > oo, tends to exp [—mot], where mo is the lowest mass in the 
channel we have selected. 

To get the best estimate for the lowest energy state it is necessary to 
optimize the choice of the operators we use. We have an infinite choice 
of lattice operators with given quantum numbers. We want to get a large 
overlap with the physical wave function by making Cq large as compared to 
the other coefficients. We will like to be able to measure a signal up to large 
distance to kill any remaining contamination originated from the excited 
states. 

The various sources of errors to be reduced in these determinations are the 
following: statistical errors, finite lattice spacing errors, finite volume errors, 
effects of higher mass states. For excellent and more technical reviews of 
error analysis in lattice calculations see for example''^'^^'*'^^'^®'^° 

We would like to get a full understanding of the large separation asymp- 
totic behavior of the heavy quark potential. The main problem in the the 
numerical calculations of the string tension a is that this quantity can be 
measured quite easily only for separations of qq of very few lattice units. We 
need the hadron, i.e. the qq system, to fit in the lattice, and the lattice spac- 
ing to be small enough for the discrete nature of the lattice to be irrelevant 
in the interaction, in order to exhibit continuum physics. But we have to go 
to very small values of in order to have a system behaving according to 
continuum physics. 

So, the situation is not easy. For small coupling constant and small lattice 
spacing we do not expect to detect the correct asymptotic behavior of the 
potential on a too small lattice. But for a too large lattice spacing we do 
not expect to be dealing with a theory close to the continuum one. To check 
scaling one has to go to high (3 values, but the lattice spacing shrinks when 
(/ — > 0. Operators of a fixed lattice extent become smaller in physical size 
and their projection over the ground state goes to zero. 

Then to get best estimate, in the continuum limit, for the lowest state 
(which for large euclidean time dominates the correlation functions) we need 
good operators. Good operators means here that we want them to have 
a large overlap with the ground state wave function. We want to be able 
to measure correlation functions up to large distance, in order to control 
the contamination of the ground state from excited states with the same 
quantum number. 

On the other hand the problem to extract the rate of an exponential fall- 
off of correlation functions of order one is a very difficult one. This is, indeed, 
one of the cases in which we can clearly see that a powerful computer is not 
really useful if we do not use the right methods to have a good control of 
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finite volume and finite euclidean effects. 
6.1. MCVM 

An alternative to concentrating on small distance correlation functions is to 
study the correlation of Wilson loops of many types, more complex than 
simple plaquettes. We will extract the true long range behavior due to the 
lowest lying glueball state. In general one needs means of modeling the 
gluon flux that propagates between the static sources, thereby reducing the 
time propagation necessary for the effective decay of higher excited states; 
in the conventional approach one is modeling this flux distribution with only 
a straight path and consequently one has to consider correlations over larger 
time separations. 

The idea is then to replace the limit of infinite euclidean time with the 
limit of infinite number of operator to pick up the ground state. In the case 
of the glueball spectrum this can be performed on different representations 
corresponding to different spin states. Then the method is to expand the 
glueball state in terms of certain Wilson loop operators acting on the vacuum 
and varying the coefficients in this expansion until the state with lowest 
energy is reached. In particular the spin glueball states are contained in 
the small plaquette operator and in infinite other large operator, of different 
shapes. One uses this basis, which may be classified according to the discrete 
group of rotations, in order to maximize the projection over the fundamental 
state. 

This procedure, named Monte Carlo Variational Method (MCVM) was 
proposed by Wilson'''', first applied for the glueball masses in ''^'''^ and then 
used in ^^^e potential function of a gluonic string with fixed 

ends in 

Let us explain Wilson's proposal applied on the lattice string states, with 
some details. 

One studies the gluonic field between a static quark and antiquark by 
considering, in principle, the traces of an infinite number of Wilson loops of 
the type shown in fig. 6. The time directed links represent static charges, 
while the space directed paths Pi and Pj create and annihilate gluons from 
the vacuum. 

The correlations Cij{R^ T) are formed by parallel transporting the path Pi 
through time T and taking its trace with the path P^. Their dependence on 
euclidean time is given by the appropriate eigenvalues in the transfer matrix 
formalism® . 

A matrix variational technique will give, for each symmetry of the gluon 
field, the linear combination of Cjj's which decreases least rapidly allowing 
to extract the largest eigenvalues of the transfer matrix from relatively small 
values of T . For example with T = a, 2a, 3a, and 3 independent path combi- 
nations it would be possible to extract 4 eigenvalues, whereas a direct study 
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Figure 6: The Wilson loop giving the correlation Cij{R^ T) between the spa- 
tial paths Pi and separated by time T as illustrated in ref.^"*. 

of rectangular loops of length up to T = 8a would be needed to obtain the 
same information in principle^^. Since loops with large values of T have small 
average values, the errors are relatively large and so accurate numerical cal- 
culations of their averages are very time consuming. Conversely the MCVM 
measures a large number of Wilson loops of different shapes (the overlaps 
between path Pi at T = and path Pi at T) at relatively small T values. In 
this way it enables a more accurate extraction of the ground state and then 
of the interquark potential and of the string tension. 

We note that if the theory has a positive definite transfer matrix the 
estimate we get at finite time for the ground state is an upper bound to the 
true asymptotic value. 

The main problem with this method is due to the statistical nature of our 
knowledge of how much we need of each different operator. The minimization 
procedure, difficult if the error on the measured operators is not very small, 
can in some cases result quite ineffective. 

6.2. The source methods 

The search of statistical improvements needed in order to get sensible mea- 
surements of long distance correlation functions has been one of the main 
developments of lattice gauge theories. One has started to deal with non 
local objects, most likely spread out over a complete time slice and a varia- 
tional calculation will need too many loops; then it is necessary to enhance 
the signal at large separations. 

The use of a source method has been proposed, in different forms, in ref. 

83,84,85,86,87 

In order to study correlation functions in a medium (the lattice filled up 
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Figure 7: The response of Polyakov loop operator, (P(t)), to a cold wall 
source, from ref. 

with gluons for example) it is possible to perturb the system and to look at 
the response to such a perturbation. This is analogous to what one does for 
studying the heat conduction properties of a given material. 

In particular one can look at the decay of the signal coming from a cold 
waW'^'^^. The source method, as applied in ref.®'', is based on the study of 
the response of the expectation value of an operator on several space slices 
at times T = 2, 3, L setting all spatial links at time equal one to the 
identity. 

Thanks to the effect of the wall which generates a strong signal one can 
measure the exponential decay for larger distances than using loop-loop cor- 
relation functions (see fig. 7). This is the reason for the use of a source in 
computations of the string tension. 

The difference from other methods is that in the case of the source the 
transfer matrix is not more positive definite and the value estimated for the 
ground state at finite time is not an upper bound to the true value anymore. 

Moreover, in the presence of the source, both the procedure of smearing 
(discussed later) the gauge field in space®'' and the variational method®^ turn 
out to be useful. 

6.3. The difference method and Langevin equation 

In the Monte Carlo approach correlation functions are obtained as expec- 
tation values of operators averaging over a suitable number of uncorrelated 
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configurations. Since large distance correlation functions are small, and they 
are normally determined by measuring terms of order one, they are dramat- 
ically affected by statistical fluctuations. For this reason a large number 
of configuration is needed to obtain a satisfactory precision in correlation 
functions measurements. 

Then the idea of perturbing the system far away from equilibrium in a 
limited space-time region and measuring the decay of correlations from such 
a zone can be pushed far away by using specific properties of the dynamical 
updating algorithm. The selection of such an algorithm can play a crucial 
role in helping to get more accurate results. 

The method based on this idea looks at differences of the values of the same 
quantity on two configurations. If we are able to keep the two configurations 
very close each other, the contribution of the statistical fiuctuations affecting 
the two configurations cancels in the difference. 

To describe the difference method^^'^^'^^'^'^ one considers two gauge systems 
K and K'{i) on two lattices of identical size X L; K is kept at inverse 
square coupling /3; K\i) has a "time slice" z = 1, 2, 3...iv) where the inverse 
square coupling takes the value fi -\- 5(3. 

We define E[j) as the average energy of the time slice j for the system 
A", and W{i^j) as the average energy of time slice j for the perturbed sys- 
tem K'{i)- Then we can show®"^, using Wilson action, that the connected 
correlation function at distance d of the energy operator C{d) is given by: 



It is easy to see that if we generate independent configurations by using 
a normal Monte Carlo simulation we are not in good conditions to compute 
correlations functions by this technique. In fact, in order to avoid noisy cor- 
relations the two sets of configurations which produce the averages K and 
K' must be similar (we are forgetting for a moment about gauge invariance). 
This cannot be achieved by a Monte Carlo algorithm, like a Metropolis up- 
dating, which does not give a continuous trajectory in the phase space; it is 
an 'yes or no' procedure and the evolution trajectories of the two systems can 
just choose if being identical or completely different. In fact a given change 
of a gauge variable can be accepted at (3 while the corresponding at (3 -\- 5(3 
is discarded. Then we need a method to generate configurations continuous 
in the /3-variable and the increase in accuracy over the Monte Carlo method 
is achieved just by allowing for a coherent cancellation of statistical errors 
between two highly correlated stochastic processes. 

Both Langevin^°'^^'®^'®^'®"^'^^ update scheme and an heat bath method made 
continuous®'' can generate sets of similar configurations. In the first scheme. 



C{d) 





5(3 



for any j such that d = \j 



— I 
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if S is the action, the gauge field configurations can be obtained as solution 
of the following Langevin equation: 



We note that two different times come into play: one is the physical euclidean 
time, one of the four dimension of the lattice, completely homogeneous to 
the other 3 space dimensions. The other, a computer time, is the evolution 
time of the unphysical dynamics of the differential equation used to formu- 
late Langevin algorithm for lattice gauge theory. After a certain amount 
of this computer time, the system will be representative of the probability 
Boltzmann distribution; therefore the Langevin equation provides a way to 
generate this distribution analytically and numerically. 

One has to be careful in noticing that the discretized version of Langevin 
equation requires the introduction of new parameters, that will have to be 
tuned in order to get good performances. 

On the other hand it is known that very efficient updating algorithms 
are the heat bath ones^"^'^^'^'*'**. It is very easy to make them to satisfy 
the requirement of a continuum trajectory; for the SU{3) gauge theory the 
method proposed in ref.** has been discussed and applied in®''. In particular 
this second method exhibits two relevant improvements with respect to the 
numerical solution of the Langevin equation: it is about ten times faster in 
reaching thermal equilibrium and does not need arbitrary parameters arising 
from the discretized version of the equation. 

Besides we note that gauge invariance effects increase the noise on the 
system and then gauge fixing or something equivalent must be imposed both 
in Langevin and heat bath updating algorithms®"^'®''. In particular the main 
disadvantage of this second method is the increasing of the noise for large 
upgrading times. An increase in the noise by a factor 10 or more occurs in 
a few iterations, then the two configurations become quite different after the 
onset of the strong noise. Moreover, using an heat bath algorithm, it is con- 
venient to measure correlation functions for operators which are a functional 
of the field smeared in space. 

6.4. The DLR and multihit methods 

Another strategy for reducing the statistical errors in the evaluation of Wilson 
or Polyakov loops for large lattices is to replace the original observable with 
its local average. One evaluates the thermal average of a link with respect 
to its nearest neighbor links. The method proposed in ref. ^^'^"^ is based on 




(6.3) 




(6.4) 
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the identity derived by Callen a particular case of the DLR equation^^. 
The idea is very simple: if two observables A and 5, with statistical errors 
respectively (A^)c and [B'^)c, are such that: 

(A) = {B) but (A^), » (5^), , (6.5) 

it is obviously much better to measure B instead of A. In ref. the SU (3) 
gauge link variable U is replaced in the Polyakov line with 

V = I dU U exp-(3H/ j dUexp-(3H , (6.6) 

where the integration is done over the links neighboring the original one. 
A heat bath method is used to evaluate U numerically. In particular 
= 14 heated hits (a multi-hit method 94,56,97-j j-^g^yg been done and it has 

been seen that the statistical fluctuations decrease with N down to a plateau 

value which is reached for N (9(10). 

6.5. Smearings and fuzziness 

The "first generation" of string tension calculations used small loops, such as 
the plaquette. It is known that for local operators the noise-to-signal ratio 
diverges as the continuum limit is approached (as discussed for example in 
ref. 98). 

In particular to check the scaling one has to compute perturbatively the 
dependence over (3 at the critical point of the theory. Then it is necessary 
to go to higher and higher (3 values despite the fact that the lattice spacing 
shrinks when — > 0. Operators of a fixed lattice extent become smaller in 
physical size, and their projection over the ground state goes to zero. To 
improve one must resort to non local operators. 

Actually several methods^"^'^^'^® are known to work to kill the unphysical 
short wavelength fiuctuations. In particular two iterative methods : the 
smearinff'^'^^^'^^^'^^'^'^^^ and the fuzzy /oops^'''^^'^^^'^™, have been introduced. 

The smearing procedure as originally proposed for SU (3) in consists in 
the construction of correlation functions for operators which are a functional 
of the field smeared in space and not in time. To this end it is possible to 
construct a block gauge field with a simple gauge invariant procedure. 

For each link of a generated configuration we consider the product of 
the other three link variables defining with it a plaquette and we sum these 
product over 4 choices which define plaquettes orthogonal to the time axis; 
the resulting matrix, projected on the gauge group, will be the new link 
variable; a smart description^^ of the procedure is in fig. 8. 

The smearing coefficient e is chosen by optimizing the performance of the 
method^°^. 

This procedure can be iterated (see fig. 9) to obtain a field more and more 
smeared in space at fixed time. One can evaluate a series of operators in 
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Figure 8: Description of Ape's algorithm for smearing from ref. 

which short distance fluctuations are more and more suppressed as the gauge 
field on the link represents a smeared average over larger and larger neigh- 
borhoods. Then one can have a set of operators, to implement a variational 
calculation, allowing to evaluate if a given euclidean time distance is asymp- 
totic; otherwise, as the amplitude of different exponential contributions de- 
pend on the operator, the ground state mass estimate from that given time 
distance can depend over the choice of the operator. 

In momentum space the effect of the smearing is equivalent to the appli- 
cation of the factor exp — A;^ (as it is shown in ref.^°^ for the example of a 
scalar field). This consideration makes more explicit the role of suppressor 
of short wavelength fiuctuations that the smearing has. 

The fuzzy procedure^^ is depicted in fig. 10 

In principle the coefficients are tunable even if ai = is chosen in 
ref. ^™ to normalize the iwLzy link back into S'[/(3). We note that the iwLzy 
procedure is inspired by the Monte Carlo renormalization group (which we 
will discuss later) methods involving factor-of-two blocking. 

Both fuzzy blocking and smearing procedures are performed only in the 
spacelike directions because unambiguous identification of the masses relies, 
through the transfer matrix, on operators with have support on one time 
slice. 

A consequence of fuzzy technique is that links paths grow exponentially 
during the iterations and there are 2^ times fewer fuzzy links than original 
links. Then, after iteration, a simple loop of fuzzy links is a complicated 
linear combination of loops of original links. Finally elementary loops of 
fuzzy links are quite non local when expressed in terms of the original links. 

A slight variation of the APF-style smearing proposing^°^ an alternative 
way of constructing operators improves the signal for a fixed number of 
smearing steps. Nevertheless, as a — > 0, all these methods become critically 
slowed down. 

An accurate test on efficiency of various smearing algorithms and a com- 
parison with the cold wall source method are performed in ref.^°^. 
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Figure 9: The string tension a estimated at different times versus the oper- 
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Figure 10: Description of Teper's algorithm for fuzzy links from ref. 
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6.6. MCRG method and Scaling violations 

The Monte Carlo Renormalization Group method is known to give accurate 
and rehable information about the critical properties of statistical systems^°'' 
108 109 110 lattice gauge theories iii'ii2,ii3_ -g ^ combination of the prin- 
ciples of Monte Carlo simulation^ with those of real space renormalization 
group ii^'iis. Here this method is analyzed only in connection to the study 
of the /3-function related to the measures of the string tension and its scaling 
violations. 

It is in general an open problem how big a lattice in a gauge theory simu- 
lation should be to do physics really near the continuum. The determination 
of the /3-function controlling the relation between the coupling constant and 
the cut-off is very relevant here. One would like for example to know where 
the /3-function changes from being representative of the non-perturbative to 
the perturbative theory. 

Since a confirmation of the approach to the continuum limit has basic 
importance in lattice QCD, systematic scaling analysis at large values of (3 
are inevitably required. In particular the measurement of the /3-function is 
crucial in lattice gauge theories. It is very difficult to establish the connec- 
tion to the continuum theory without it. The /3-function describes how the 
bare coupling should be tuned in order to keep all the physical predictions 
independent of the cut-off in the continuum limit. The /3-function is unique 
in this sense but not universal: it is different in different lattice formulation. 
In particular the /3-function depends on the lattice action chosen. Only the 
two leading terms in its perturbative expansion are universal: 

m = -Po9'-f3ig' + 0{g') . (6.7) 

as given in eq. (4.5) For large cut-offs (small bare coupling values) these terms 
dominate and define a universal scaling behavior, the asymptotic scaling. 
Outside this region, but still in the continuum limit, the scaling behavior 
is described by the full, and in general unknown /3-function. This function 
gives the way asymptotic scaling is approached, connects numerical studies 
with perturbation theory, reveals the existence of possible phase transitions. 

In the MCRG approach, not the /3-function itself , but a related quantity 
A/3 = A/3(/3) is determined, which gives the change of the coupling 

f3 ^ f3 - A(3{(3) . (6.8) 

when the (dimensionless) correlation length (or the cut-off) is decreased by 
a factor of b. Here b is the basic change of scale in a single renormaliza- 
tion group (RG) step (6 = 2 in the following). At the couplings (3 and 
fil = fi — A(3[(3) the model has identical long-distance properties, only the 
(dimensionless) correlation length ^ differs by a factor of 2. 
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A central part of any MCGR study is the block transformation (BT) 
which calculates the block fields on the next coarser block lattice from av- 
eraging link variables over appropriate regions of the current lattice. The 
system under consideration is divided into blocks and a smaller number of 
block variables are defined by averaging in some fashion over the original site 
variables. By studying the way in which the site hamiltonian (or lagrangian 
for lattice gauge theories) "fiows" into the block renormalized hamiltonian 
one can determine the critical properties of the system. 

The most important properties of a BT are: 

1. The partition function of the two systems (unblocked and blocked) is 
unchanged but the correlation length of the block variables in units 
of the block lattice spacing is one-half that of the original variables in 
units of the original lattice spacing. This rescaling of lengths is the key 
to any RG transformation. 

2. The interaction between the blocks is described by a new action which 
in general will contain all kinds of interactions. 

Here we have used the crucial assumption that a fixed point exists and that 
it is short ranged. Renormalization generates an infinite number of couplings, 
but we assume that a small number of short range points is enough to describe 
the system at a given scale and the long distance physics is preserved. 

The new action generated through iteration of the RG transformation 
can be considerated as a point in the multidimensional space of the different 
couplings constants and a BT transformation is expected to have a fixed 
point somewhere in the (3 = oo hyperplane of this multidimensional space. 
A single renormalized trajectory (RT) starts form this point (see fig. 11 from 
refill^). 

Both the positions of the fixed point and of the RT are not universal and 
depend on the BT. We start with the standard action at a given (3 value (for 
large (3). The effective actions obtained after a few BT will move close to the 
RT. The same will happen if we start at some other coupling values /?': by 
tuning (3' the points of this second sequence, which lie on the RT, can reach 
the corresponding points of the first sequence but one step behind. 

Different procedures have been proposed to implement the MCRG. The 
first one originally proposed by Wilson^^^ is to find an action close to the RT 
of a given BT. A different strategy is to look for an improved BT whose fixed 
point and RT lie close to the standard action^^'''^'''^^^. 

An other method^^^ is based on the Schwinger-Dyson equations, providing 
the block renormalized couplings and the RT for any value of the initial 
couplings, also far from the critical values. An accurate review of these 
methods is in ref.^^"*. 
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Figure 11: The RT represented in the hyperplane of the multidimensional 
coupling constant space. The fixed point lies in the (3 = oo plane and the 
standard action corresponds to the points of the C2 = C3 = ... = axes 

In this section we will describe a generally applicable renormalization 
transformation (BT) for analyzing lattice gauge theories as presented origi- 
nally in ref.^^"^. Others BT are used in the literature: in fact many different 
BT can be used to analyze a given model and their comparison is useful to 
check block transformation universality^^'*'^^^ and efficiency. 

The use of real-space RG in lattice gauge theories has proven more difficult 
than for the corresponding spin systems because of the high symmetry of a 
gauge theory. Most simple definitions of block spins fail to preserve the gauge 
symmetry and can only be used after some form of gauge fixing. 

The transformation introduced in^^"^ avoids gauge fixing either globally 
or within specified blocks^^^, is directly applicable to arbitrary lattice gauge 
theories and the renormalized Hamiltonians (actions) retain the full gauge 
symmetry of the original model. 

This BT is illustrated^^"^, for simplicity, on a two-dimensional lattice in 



The lines represent the Uij gauge variables and the intersections are the 
sites of the original lattice. For a transformation with scale factor 2 only 
the circled sites remain in the new renormalized lattice. The renormalized 
gauge fields connecting two sites of the blocked lattice are constructed from 
the operator products along the paths marked A, 5, and C as shown in fig. 12 

To complete the description of scaling analysis by MCRG methods, per- 
formed by several gj-Qups^^'''^'''^^^'^^"^'^^'*'^^^'^^^ we recall two methods for cal- 



fig.l2. 
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Figure 12: Diagram of the renormalization transformation for lattice gauge 
theories as proposed in ref. sketched in two dimensions. Operator products 
are defined along each of the paths marked A, B, and C. 

culating the non perturbative /3-function directly: the loop ratio method^^®'^^ 
117,97,120 ^j^g Wilson's two lattices methodiii'ii2-ii3,i20,i23,i26_ ^j^g second 

method, designed for models with one relevant direction in coupling param- 
eter space, is used in the one of the most recent scaling study^^^ of quenched 
pure gauge lattice QCD. 

This method is based on the blocking which is repeatedly performed down 
to 2'^ for two lattices, one of size L (at fi) and the other of size iv/2 (at (3 — A(3). 
To match long range physical contents on both lattices, a set of Wilson loops 
on one blocked lattice is compared with the corresponding one on the other 
blocked lattice. 

The expectation values from the two simulations are then compared on the 
same size lattices, i.e. the ones from the larger starting lattice L blocked one 
more time than those from the smaller lattice. The test for convergence of 
the two systems is that the expectations values should match simultaneously 
at the last few levels. For early matching between blocking trajectories, the 
blocking transformation is controlled by a parameter q which governs the size 
of Gaussian fiuctuations around the maximal SU{3) projection of the block 
link variables^^"^. 

The coupling shift A/3 is determined at the value of q where the mismatch 
of the two sets of Wilson loops is minimum; planar 1x1 and 1x2 Wilson 
loops, non-planar 6-link ("twist" and "chair") and 8-link ("sofa") loops are 
measured on blocked lattices (fig. 13) 

In the frame of the MCRGM different updating methods have been used: 
the heat-bath^'* and the over-relaxed pseudo heat-bath algorithm ^^"^ . 
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7. The use of different actions 

An important question raised just after the first measurements of the string 
tension was how close were the measured values to the real, continuum con- 
tent of the model considered. 

Wilson's point of view^ that confinement persists for arbitrary small values 
of the bare charge seemed to be confirmed by the first results. A possible 
useful internal check was to compare results obtained by using different lattice 
actions in the scheme of the universality. A possibility to reduce lattice 
artifacts using different actions is based on the idea that for a large class 
of actions the physical predictions are universal while the cut-off effects and 
in general the lattice artifacts are not. No matter what the form of the 
lattice action is once the gauge group is fixed the continuum physics should 
be defined. In particular the string tension as measured by Monte Carlo 
methods should behave for small enough bare charge in a universal way, up 
to a multiplicative constant determined by the action in question. The change 
of the lattice action corresponds to modify the scheme of regularization but 
the continuum limit should not be affected. 

The corresponding (A) scale parameters can be exactly connected by a 
perturbative calculation (in the (/ — > limit) i4i, 135,147,136 -g pQggi]3}g 

to check if the Monte Carlo simulations follow such predictions. 

The results I37,i35,78,i38 153,139,157,140,154,106 referring to different actions, (for 
example the Manton's action^^®, the heat kernel action 129,130,131,132 ^j^^ 
mixed actions consisting of terms belonging to fundamental, adjoint and 
others representations^"^"^'^°^), show that the scaling form of the string tension 
is reached in different coupling region for the various actions. In ref.^"^"^ it 
is discussed how a modification of the Wilson action can introduce a first- 
order phase transition between the strong and weak-coupling regime in SU (2) 
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lattice gauge theory; nevertheless this transition is not deconfining because 
it can be continued around in a larger coupling space. The hypothesis that 
more structure could be present in the string tension was investigated also 

A test of the universality and scaling for the measure of the string tension 
whit Manton's, Symanzik's tree-level (discussed later) and Wilson's action 
in SU{2y^'^ shows that significant violations of scaling are similar for the 
three actions and suggests that the reason for this violations is common to a 
whole subclass of lattice discretizations and could be associated with the un- 
derlying cubic lattice substructure. Manton's and Symanzik's actions violate 
asymptotic scaling in the same direction as Wilson action but significantly 
more weakly; nevertheless such an improvement is payed with the difficulties 
arising from the violation of positivity characterizing improved actions^'*^. 

So the use of the improved actions had a limited effect on actual calcu- 
lations until recently. This has also been because of the poor performance 
of bare perturbation theory. A description of different methods of improving 
the action, with their advantages and disadvantages can be found in ref.^^"* 
and the state of art of these procedures in^'*"^'^'*'* and references therein. 



7.1. Symanzik's improved actions 

In the frame of finding an improved action in lattice gauge theories to re- 
duce the effect of operators that lead to scaling violations the Symanzik's 
improvement program^'*^'^'*^'^'*'''^'*®'^'*^'^^^ assumes a particular relevance. 

This program is based on building a systematic reduction of the cutoff 
dependence of physical quantities, diminishing the corrections to continuum 
theory stemming from finite lattice spacing. If a theory is regularized by 
introducing a lattice spacing a, when a — > the regularized theory should 
go to the continuum theory; however, when a is different from zero the two 
theories are different. Symanzik's improvement program consists in choosing 
the action on the lattice in such a way that the difference between the two 
theories is as small as possible. 

Symanzik's scheme specifies in a precise way how to quantify the difference 
between the lattice regularized theory and the continuum theory and shows 
that the improved action, defined following the criterion he has proposed, 
can be computed in perturbation theory or measured in conceptually simple 
Monte Carlo experiments. 

As noticed in section (4) any of the physical quantities calculated in a 
lattice regularized version of the gauge theory with perfect scaling will obey 
the renormalization group equation (4.1). 

For a generic action this will however not possible and, according to per- 
turbation theory, the best one can achieve is approximate scaling near the 
continuum limit, i.e. for a — > 0: 
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d d 

- a— — h /?((/)——) [physical quantity] = (9[a^(lna)] (^-l) 
da dg' 

The suggestion of Symanzik's improvement program is to reduce the scal- 
ing violations to terms at least of order (9(a'*(ln a) to all orders of perturbation 
theory including in the lattice action suitable chosen irrelevant terms. Actu- 
ally it is important to note that not only spectral quantities can be improved, 
but also the (properly normalized) n-point functions of the fundamental fields 
in momentum space. 

In ref.^^^ on can find a review of the theoretical basis, the possible appli- 
cations of Symanzik's improvement program and a discussion of his possible 
extensions to be compatible with the positivity and symmetry of the trans- 
fer matrix. In fact the Osterwalder-Schrader positivity condition is crucial 
in numerical simulations: it is at the heart of the MCV method discussed in 
section (6.1). 

For 4c/ SU{N) lattice gauge theories the Symanzik's tree-level improved 
action (TIA) has been determined^'*^'^'*'''^'*®'^'*^. A motivated ansatz^'*'' for the 
improved action includes Wilson loops up to length 6. The result^'*^'^'*® is a 
function also of the planar rectangular double plaquettes of size 1x2. 

Different Monte Carlo studies of string tension have been carried out using 
Symanzik improved actions^^^'^^"^'^^'*'^'*^. In particular we refer to studies of 
Creutz ratios up to x(4,4)^^'*: an improvement of previous string tension 
estimates ^"^'^^ is allowed and it is also shown that Creutz ratios xi^^ J) 
not stable under increasing {I, J) from (3,3) to (4,4). It is this a clear 
tendency of the estimated string tension to decrease as already noted in that 
period by a potential analysis^^'' and a Polyakov loop measurement^^^'^^^. 

We also note the relevance of a consistent string tension measurements 
involving both the improved action and corresponding improved Creutz ratios 
which expression at tree level order is^'*®: 

X'{LJ) = E CL^n HW{I + m, J + n)] . (7.2) 

where the improved „ coefficients and the corresponding ones of the stan- 
dard Creutz ratios can be found in ref.^^"*. 

Scaling of the improved string tension sets on at a smaller correlation 
length than for the standard plaquette action but scaling windows, for similar 
lattice sizes, look similar and therefore, from a practical point of view, the 
improvement is modest. An important consistency check is, of course, the 
observation that the improved definitions^'*® of Creutz ratio work well. 

In a following study*"*^ calculations for larger lattices show that the string 
tension follows 2-loop perturbative theory more closely in the case of the 
TIA and Manton's actions than in the case of the standard plaquette action; 
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still there is a warning on the difficulties caused by the violation of positivity 
of improved actions. In particular in ref. ^"^^ it is showed that in the case 
of the Wilson block-spin-improved action^^^ the violations of positivity make 
impossible to extract accurate masses at small (3 values. One has to be careful 
in giving a final word about the program based on understanding continuum 
physics at smaller (3 values using improved actions. 

8. String picture 

An alternative way to parameterize the lattice data in the large r region 
is to assume a simple long distance picture of QCD, the one of a chromo- 
electric fiux confined to a tube. The color fiux joining a pair of quarks in the 
confining phase is concentrated inside this tube of small but finite thickness. 
The possibility of describing the long-distance dynamics of gauge theories in 
the confining phase by an effective string theory^^®'^^^'^^° is based on the very 
intuitive assumption that the color flux connecting a pair of distant quarks is 
squeezed, in the conflning phase, within a narrow flux tube (string) carrying 
constant energy density; as a consequence, the energy of the system increases 
with the separation of the quarks. 

It is generally believed that this thin tube behaves like a vibrating string 
when the quarks are pulled very far apart. This is also supported by the 
strong coupling expansion of the lattice gauge theories which can be formu- 
lated as a sum of weighted random surfaces with quark lines as boundary. The 
string picture of conflnement in the strong coupling approximation within the 
framework of the Hamiltonian formulation was discussed originally by Kogut 
and Susskind^^^. 

The fluctuating string can be taken as the description of the color flux 
between colored sources: since the string model appears to be reasonable 
it is interesting to quantify the nature of this string. One can study the 
modes of a scalar gaussian string^^^ to derive which terms contribute to the 
potential. In particular, in lattice calculations, static color sources can be 
used, so that an appropriate method to explore the string model is a study 
of the long range static potential V{R). 

Nevertheless the action describing the effective string in the continuum 
limit is substantially unknown. The simplest possible assumption is that the 
action of the effective string is described by the Nambu-Goto string^^^'^^"^ in 
terms of [d — 2) free bosonic flelds associated to the transverse displacements 
of the string. The leading behavior of a R X T Wilson loop is^^"*: 



rji If? 1 OO 

-lnlF(i?,r) = ai?r+Ki?+r) + c-(J-2){^+^ + -^ln(l-e^)}. 

n=l 

(8.1) 

The last term is universal and only depends on the number of transverse 
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dimensions [d — 2). The parameters cr, c depend on g. 

Different Icinds of difficulties are present in applying the bosonic string 
outside the critical space-time dimension of 26. Modifications^^^'^^^ intro- 
duced to transform the Nambu-Goto action into a consistent theory make 
it very difficult to evaluate physical observables which should be the final 
goal of the effective string picture of the gauge theories. In effect, consid- 
ering only the infrared behavior of the interquark potential, the effective 
string is asymptotically described by a two dimensional conformal field the- 
ory formed by [d — 2) massless free bosonic fields which is the infrared limit 
of the Nambu-Goto action called the free bosonic string. 

Different schemes are proposed to overcome the too drastic approxima- 
tions given by the free bosonic string. A simple modification of this pic- 
ture^^'''^^® consists in a suitable compactification of the boson fields. It fulfills 
the constraints dictated by the gauge system^^^ which cannot be obeyed by 
the free bosonic string and fits well the numerical simulations of the gauge 
systems in three and four space-time dimensions and with various gauge 
groups. The effective string picture arising in this way coincides with that 
proposed previously^ in order to fit accurately the numerical data on the 
expectation value of the Wilson loops for various gauge systems in three and 
four space-time dimensions and based on a different argument^^^. 

A numerical exploration of the string picture studying the long range 
static potentials between static fundamental color sources corresponding to 
the ground state and excited states has been performed in ref.^''. The predic- 
tion of this picture is that the ground state potential between static quarks 
consists of a term linear in the separation R of the sources and of correc- 
tions to this term originating from the zero point motion of the quantized 
string. The leading correction is proportional to 1/R and the constant of 
proportionality has the universal value 7r/12^^"^ for a general class of bosonic 
string theories. Considering the Nambu string^^^'^^"^ with fixed ends in four 
dimensions the corresponding potential for the different string states is given 
by the formula: 

Vn = [K^R^ - ttK/Q + 2nNKf'\ = 0, 1, 2, 3, ... (8.2) 

The agreement between predictions of this picture and lattice simulations is 
good at relatively large values of i?^'' where the width of the flux tube is 
smaller than its length and the string model is expected to apply. Problems 
and discrepancies arise for small R. 

However the quantum nature of the string is only revealed by study of the 
string fluctuation term and more subtle tests involve the excitations modes 
of the string which have energy vr/i? per mode at large R. Then the string 
picture is clearly a useful guide to the excitation energies at large R. 

As it is well known^^"^ the string fluctuation component in V{R) for a 
bosonic string behaves as 7r/(12i?) at large R and a keen test concerns just 
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this self-energy term which is difficult to separate unambiguously from the 
contribution to the potential of an effective Coulomb term dominated by one- 
gluon exchange also behaving as 1/ R. Suggestions have been made that other 
self-energy expressions are appropriate for different string models. The Dirac 
string described above^^'''^^® (involving the Coulomb piece fixed by lowest 
order perturbation theory) leads to a self-energy term four times smaller. 

So it would be useful to determine the self-energy contribution accu- 
rately from lattice calculations. The difficulty to separate this term from 
the Coulombian component lead to the proposaP''^ to explore 3-dimensional 
gauge theories where the Coulomb and string fiuctuation effects are clearly 
different [logR and 1/R respectively). Attempts to separate effects due to 
gluon exchange (Coulomb) and string fluctuations (conflning force) from an 
analysis of the spin orbit potentials'''"^ have been performed but more 
accurate results are needed to conflrm the validity of this approach. 

Another way is to study closed strings via states of electric flux which 
encircle the periodic spatial boundary conditions, called torelons, of size L. 
These are systems which energy E{L) will again have a string tension com- 
ponent KL, but where the string fluctuation term behaving as tt/[3L) for a 
bosonic string is unaffected by any Coulomb component; so a search for the 
self-energy is propitious. To extract the string fluctuation component one 
should consider long strings so that the string model is most likely to apply; 
L X X T lattices with both S and T large have been analyzed'''^. 

Accurate measurements of the energy of the ground state torelon E[L) 
show conflnement generated through a color flux tube having energy which 
increases with his length: the existence of a constant conflning force (the 
string tension K ) is conflrmed in a non-abelian SU (2) color gauge theory. 
A careful study^''^ of E{L) of the torelons can expose the possible self-energy 
contributions in a hadronic string. Following ref.^'' the parameterization is: 

Eiiy = {KLf - ^ , (8.3) 
which for large L becomes: 

E{L)/L = K - fir/SL'' . (8.4) 

The comparison'''^ of the flts with the data is illustrated in flg.l4. 

The slope in the flgure is directly proportional to /; the continuous line 
shows the flt which gives strong evidence for a term in E[L) with behav- 
ior at large L. This term is expected from a string self-energy contribution 
and the coefficient /, compared with that of simple string models, is found 
completely consistent with the bosonic string value of / = 1 and inconsistent 
with / = 1/4, the Dirac string value. For the conventional string model there 
is no discrepancy between the string tension extracted from potentials and 
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Figure 14: The torelon energy E[L) versus L for a periodic boundary of size 
L. 

torelons. This results are valid for pure gauge SU{2) color fields and similar 
conclusions can be expected for SU{3). 

In conclusion the hadronic string model is a useful guide to the nature of 
the extended color fiux but evidence for the characterization of self-energy 
from a particular type of string is not yet completely compelling. 

9. Conclusions 

We have discussed the status of the results of numerical measurements of the 
string tension in lattice gauge theories, in particular for SU{2) and SU{3). 
This relevant physical quantity has a good scaling behavior^^^'^°. In order 
to make that explicit we have to use an improved coupling constant. In fact 
the string tension remains constant, for example in SU{3)^ in the (3 region 
5.70 — 6.80 where the lattice spacing changes by a factor 6. Although the 
quality of the agreement varies from one frame to the other, and a 20 percent 
discrepancy between different methods^^'^°'^"^ remains, the approach to the 
expected continuum limit seems confirmed. 

Present and future efforts are and will be directed to the study of larger 
lattices, at larger values of the coupling constant. The other main direction 
points to a better study of full QCD. 

In fact the confinement in pure gauge QCD is visualized as the formation 
of a color electric fiux-tube between two isolated static quarks. The strength 
of the linear contribution is connected to the string tension a. This picture 
represented by the phenomenological spin-independent potential as given by 
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the eq.5.3 changes when including dynamical quarks. At short distances dy- 
namical quarks renormalize the gauge coupling changing the value of aqq[r). 
At long distances the creation of a q — q pair from the vacuum becomes 
favorable, and the string breaks. 

This review has been discussing only quenched lattice gauge theories, but 
we note that lattice calculations that include some sea quark effects begin 
to exist (see for example ref. I82,i80,i8i references therein); moreover 
a program to estimate these effects phenomenologically from quarkonium 
spectroscopy^'''' has already given promising results for a full QCD determi- 
nation of the strong coupling constant ^''^ including an even more realistic 
spectrum of sea quarks. 

In this new context all the efforts to establish effective algorithms have to 
be replicated. One has to do a new tuning and to try and select the most 
efficient algorithms. The high (3 limit will involve a fight against critical 
slowing down and small quark masses call for efficient methods to inverse 
large sparse matrices. For example in we have discussed a very promising 
approach to the measure of properties related to confinement in fully coupled 
theories. We note that also the method determining directly the strong 
coupling from the interquark force^° is applicable to full QCD. 
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